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We show how discrete squeezed states in an iV 2 -dimensional phase space can be properly constructed 
out of the finite-dimensional context. Such discrete extensions are then applied to the framework 
of quantum tomography and quantum information theory with the aim of establishing an initial 
study on the interference effects between discrete variables in a finite phase-space. Moreover, the 
interpretation of the squeezing effects is seen to be direct in the present approach, and has some 
potential applications in different branches of physics. 
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^ ; I. INTRODUCTION 

The first prototypes of continuous wavepackets related to the coherent and squeezed states appeared in the literature 
i in 1926 through the publication of three seminal papers by Schrodinger [l[, Kennard Q, and Darwin 0], where the 
' evolutions of Gaussian wavepackets for an oscillator, a free particle and a particle moving in uniform constant electric 
, and magnetic fields were treated for the first time. In fact, they represent not only an initial historic mark in the study 
of nonclassical states of the electromagnetic field, but also, from the contemporary point of view, the cornerstones of 
modern quantum optics p|. Since then a huge number of works dedicated to application, generation, and detection 
of nonclassical states have been published in many specialized scientific journals covering different areas of knowledge 
in physics (such as, for instance, solid state physics, nuclear physics, high energy physics, general relativity, and 
cosmology). Nowadays, beyond these fundamental features, these states have also potential applications in quantum 
information theory and quantum computation since continuous- variable entanglement can be efficiently produced using 
squeezed light and nonlinear optics [falfj- However, if one considers physical systems with a finite-dimensional space of 
states, a sound theoretical framework must be employed in order to describe properly the nonclassical states. In this 
sense, it is worth noticing that there are various formalisms proposed in the literature for finite-dimensional Hilbert 
spaces with convenient inherent mathematical properties which can be applied in the description of these quantum 
systems @, H, H, floL [Til EE EE EE EE EE EE HE EE]- In particular, let us focus our attention upon the approach 
developed in Ref. [11] for the discrete representatives of the coherent states which has it s alg ebraic structure based on 
the technique of constructing unitary operator bases initially formulated by Schwinger [20] - Basically, this approach 
consists in proposing a new specific operator basis (whose discrete labels are congruent modulo N) which allows the 
direct mapping of Weyl-Wigner representatives in an A^ 2 -dimensional phase space through a simple trace operation; in 
addition, a generator of discrete displacements for this particular phase space, which exhibits an additional topological 
phase related to the mod(iV) invariance, is also derived and its properties properly discussed in details. Subsequently, 
we presented in Ref. [181 ] an ab initio construction that inherently embodies the discrete analogues of the desired 
properties of the Cahill-Glauber formalism [2l| . where some direct applications in quantum information processing, 
quantum tomography, and quantum teleportation were explored. Furthermore, this specific construction process 
opens the possibility of dealing with discrete wavepackets whose widths are modified, in principle, by parameters that 
mimic the continuous squeezing effect. 

The main goal of this paper is to fill this breach by constructing a consistent formalism for the discrete squeezed 
states in finite-dimensional spaces. For this purpose, we first present briefly some useful results obtained in Ref. [H| 
for finite-dimensional Hilbert spaces with special emphasis on the mod(./V)-invariant operator basis and the discrete 
counterpart of the nondiagonal projector related to the coherent states. Next, we define a wavepacket for the discrete 
squeezed- vacuum state in connection with the equivalence relation between the wavepacket phase space representation 
and the phase space generated by the continuous squeezed states [22| . Consequently, this definition leads us to study 
two special classes of discrete squeezed states in an ^V 2 -dimensional phase space and establish some results for their 
respective diagonal projectors in complete analogy to the continuous ones. Such discrete extensions are then applied 
to the context of quantum tomography and quantum information theory with the aim of enlarging our knowledge on 
the quantum interference effects in a discrete phase space. In particular, we attain new results within which some of 
them deserve to be mentioned: (i) we show that a generalized version of the scattering circuit introduced in Ref. [23j 
can be used to measure any discrete Wigner function and/or its corresponding characteristic function in the presence 
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of squeezing effects; (ii) we establish a first reasonable measure of functional correlation between the discrete variables 
of an iV 2 -dimensional phase space which can be extended within the context of quantum information theory, as a 
measure of entanglement in multipartite systems (24], or in other branches of physics whose systems of interest can 
be described by finite state spaces; and finally, (iii) we apply this measure for the discrete squeezed-vacuum state and 
verify that its maximum value depends essentially on the Hilbert-space dimension. 

This paper is organized as follows. In Section II we discuss briefly the extension of the Cahill-Glauber formalism 
for finite-dimensional Hilbert spaces with emphasis on some essential features exhibited by the mod(iV)-invariant 
operator basis. In Section III we introduce the discrete representatives of the displaced squeezed- vacuum states and 
squeezed coherent states in a finite phase-space, where their respective diagonal projectors are derived in details. The 
results obtained are then applied in Section IV within the context of quantum tomography and quantum information 
theory in order to establish two different but complementary methods for studying quantum interference effects on 
such phase space. Finally, Section V contains our summary and conclusions. Furthermore, we add two mathematical 
appendices related to the calculational details of some basic expressions used in the previous sections. 



II. PRELIMINARIES 



The cornerstone of the extended Cahill-Glauber formalism for finite-dimensional spaces is given by the mod(iV)- 
invariant operator basis |I8] 
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which consists of a discrete double Fourier transform of the extended mapping kernel S( s \r),£) = \Jtf{r\ 1 £)] -s S(7y, £), 
where the labels r\ and £ represent the dual momentum and coordinate variables of a discrete iV 2 -dimensional phase 
space. In particular, these labels are congruent modulo N and assume integer values in the symmetrical interval 
[—£,£] for I = (JV — l)/2. The extra term J(r(r],£) is defined in this context through the ratio ./#(?7, £)/^#(0, 0), 

= ^[^3(a?7|ia)^3(ae|ia)+^3(a?7|ia)^4(ae|ia)e im ' 

+ i?4(ar7|ia)tf 3 (a£|ia)e i7r « + ^ 4 (o77|ia)^ 4 (ae|ia)e i ' r( " + « +JV) 

being the function responsible for the sum of products of Jacobi theta functions evaluated at integer arguments [2|| , 
with a = (2iV) _1 fixed. Furthermore, the complex parameter s obeys |s| < 1, and 
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corresponds to the symmetrized version of the unitary operator basis proposed by Schwinger [20j | . Finally, it is worth 
mentioning that a comprehensive and useful compilation of results and properties of the unitary operators U and V 
can be promptly found in Ref. [ll| ; here the main focus of our attention is related to the essential features exhibited 
by Eq. CD. 

For instance, the set of N 2 operators {T'*'(/i, v)~\ n,v=-i,...,l allows us to decompose any linear operator O through 
the expansion 



= ^ fl ( - s W)T (s W), 



(2) 



where the coefficients 0( _s ) (/x, v) := Tr[T( -s ) (/j,, v)0] are associated with a one-to-one mapping between operators 
and functions belonging to a finite phase-space characterized by the discrete labels [i and v. The first practical 
application of this decomposition is related to the calculation of the mean value 
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which leads us to define the parametrized function F^(/i, v) :— Tr[T( s )(/i, i>)p] as being the double Fourier transform 
of the discrete s-ordered characteristic function E^ s \r], £) := Tr[S( s )(n,£)p], namely 



exp 
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Thus, in complete analogy to the continuous case, the discrete Husimi, Wigner, and Glauber-Sudarshan functions can 
be directly obtained from Eq. Q for specific values of the compl ex p arameter s. For physical applications associated 
with quantum tomography and quantum teleportation, see Ref. [la ]. 

After this condensed review, we will obtain a decomposition for the nondiagonal projector related to the set of 
discrete coherent states v)}u,u=-t,... £■> whose mathematical properties were established in Ref. For this task, 
let us initially consider the expansion formula for O = v){n', v'\ and s = — 1 fixed. The next step consists in 
substituting the mod(iV)-invariant operator basis T'- -1 ^/!, v) into Eq. ((2|) and evaluating the sums over the discrete 
indices p, and v. After lengthy calculations, the analytical expression for the nondiagonal projector assumes the exact 
form 
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being reduced to the diagonal projector T^ 1 ) (//, i/) = u){fj,, v\ when [i! — \i and v' = v. Consequently, the 
probability of finding the discrete coherent state v) in the nth excited state - described in [2(|[27| by a complete 
set of number states {|n)} n =o,...,JV-i - can be promptly estimated as follows: 
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where J^(n, £) :— „#„(n,£)/./#„(0, 0) is connected with the diagonal matrix element (n|S(?7, £,)\n) through the equality 
Jfr n (v,0=VN{n\S{r 1 ,0\n}, l md 



^ n (r),0 :=V^ — (l + 23)- 1 M(7 7 ,^; 3 ) 
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In this definition, Ji n (j},K) depends explicitly on the n-th order derivative of the product (l + 2j) 1 M(?y,^;3) evaluated 
at the point 3 = 0, 

+ ^4[a7 7 |/(3)]^KI/(3)]e i ^+^N|/(3)]^[ae|/(3)]e iT( '' + « +JV) } 

being an auxiliary function whose analytical expression differs essentially from ^{rj^) by the presence of /(3) = 
ia(l — 23) (1 + 23 ) -1 in the second argument of the Jacobi theta functions (see appendix A for technical details). Note 
that for n — (reference state) the function J£o(?7,£;) coincides exactly with J^(r],£) := (r),£\Q,0), as it should be in 
this case. Furthermore, adopting the mathematical procedure established in Ref. [28[ for the continuum limit, we 
also verify that P n (/x, v) goes to 



P n (q,p) = -7 exp 



q 2 + p 2 
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in the limit N — ► 00. In the next section, we will focus upon the following question: 'How can we propose a family 
of squeezed states in a finite-dimensional Hilbert space by means of a discrete wavefunction whose width is modified 
through a squeezing parameter?' 



III. DISCRETE SQUEEZED STATES FOR FINITE-DIMENSIONAL SPACES 

The reference (vacuum) state |0) has an important role within the theory of continuous coherent states [29j since 
they can be generated by means of a unitary displacement operator T>(q,p) acting on this reference state, namely 
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\q,p) ■— T>(q,p)\0). Pursuing this line, we have shown in Ref. [11| that the Schwinger operator basis elements 
{S(?7, £)}?7, a l so ac ^ as displacement operators on a particular reference state to generate discrete coherent 
states, according to the definition \fi, v) :— ^/NS(v, — /i)|0). In this formalism, the normalized discrete wavefunction 
associated with the reference state in the coordinate-like representation {|it 7 )} 7= _^.. , t t is expressed as 



(«rlO) = 



2q 



.#(0, 0) 



1/2 



T?3(2o 7 |2io), 



(7) 



where the width of the Jacobi ^-function has a constant value equal to v2a. Now, let us suppose that there exists 
an analogous wavepacket for the discrete squeezed- vacuum state |0;s) whose formal expression reads as [30j 
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with ^# s (0, 0) given by 
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In principle, the squeezing parameter s is positive definite and assumes any values in the interval s € (0, oo). Indeed, 
the introduction of this factor in the second argument of the "^-function allows us not only to modify the width of the 
wavepacket ([7]) for values greater or smaller than \/2a, but also to generalize the theory of discrete representations for 
finite-dimensional Hilbert spaces in order to include the discrete squeezed states within the scope of a general theory 
from which the coherent states represent a particular case. Moreover, note that (u 7 |0;s) is a discrete representative 
of that equivalence relation established in Ref. [22j between the wavepacket phase space representation and the phase 
space generated by the continuous squeezed states, and represents the cornerstone of the present approach. Next, we 
will study two special classes of discrete squeezed states in an A^ 2 -dimensional phase space following the mathematical 
recipe described in the previous sections. 

A. Discrete representatives of the displaced squeezed-vacuum states 

Let us initially decompose the projector |0;s)(0;s| through the expansion formula for s = — 1 fixed. After some 
algebraic manipulations, it is easy to show that 



|0;s)(0;s| 
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has a similar mathematical structure to that obtained for the projector of discrete coherent states - here represented 
by the operator T^ _1 ^(0,0) - where, in particular, J^(ry,£) := ^ 3 (r], £)/„# s (0, 0) can be promptly determined with 
the help of the auxiliary function 
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The next step consists in applying on both sides of equation © the generator \/NS(v, — /x) with the aim of obtaining 
the discrete representatives of the displaced squeezed-vacuum states. Such procedure leads us to establish a closed- 
form expression for \fi, i/;s)(fi, u;s\ by means of the identity [Tlj 
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FIG. 1: Three-dimensional plots of (a,b,c) Po(fi,i / ;s) and (d,e,f) Pi(fi,v;s) versus {/j,, v} € [—8,8] for N = 17 fixed, and 
different values of squeezing parameter: (a,d) s — 1 (discrete coherent states), (b,e) s 2 — 5 and (c,f) s~ 2 = 5. Note that both 
the squeezing and stretching effects occur on the 17 2 - dimensional phase space labelled by the dimensionless discrete variables /j, 
and v. Indeed, these effects are directly linked to the modified second argument of the Jacobi theta functions, which corroborate 
the approach here adopted for the discrete version of the displaced squeezed-vacuum states. 



Hence, Eq. (|10p can be considered as being the first analytic discrete expression of the displaced squeezed-vacuum 
states in a finite phase-space. 

An alternative form of Eq. (fit))) can also be found by expressing the product N~ 1 ^ 2 J(f{ri, £)S(?7, m terms of a 
double Fourier transform of the diagonal projector v'\ for the discrete coherent states, that is, 
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where K(?7,£;s) := J(f s (r],£)/J&'(T},I~) denotes a new kernel-like function which is responsible by the squeezing propa- 
gation into a discrete finite-dimensional phase-space. An immediate application of this equation refers to the evaluation 
of \{n\n, v\s)\ 2 through the results obtained until the present moment. For such task, it is sufficient to calculate the 
overlap probability distribution P n (fx, v;s) by means of the equation 
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which describes how the squeezing effect is propagated for a given initial distribution function P ra (/i', v 1 ). Next, we 
substitute Eq. ^ into this expression and carry out the sums over //' and v 1 in order to obtain the compact form 
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N 



exp 



»?,«=- 



27ri / , ^ 



^,0^,0- 



(12) 



6 



This result states that P n (fi,v;s) is promptly estimated through the double Fourier transform of the product 
J^siVi 0^n( 7 ?; £)> an d for fi = 1 the overlap probability P n (n,v) can be immediately recovered. 

Figure 1 shows the three-dimensional plots of P n (/J,, v; s) versus {fi, v] E [— £,£] for (a,b,c) n — (vacuum state) 
and (d,e,f) n = 1 (first excited state) with N — 17 fixed. In the numerical investigations, we have adopted some 
specific values for s in order to illustrate the squeezing and stretching effects within a finite-dimensional phase space 
characterized by the discrete variables \x and v. For instance, pictures (a,d) correspond to the particular value s = 1 
(absence of squeezing effect) and represent the three-dimensional plots of P n (/i, v) - see Eq. ^ for technical details 
- when n = and 1, respectively. On the other hand, pictures (b,e) are related to the value s 2 — 5 and show some 
typical effects that confirm the quantum assignment of the displaced squeezed-vacuum states. Furthermore, pictures 
(c,f) also exhibit similar structures for s~ 2 = 5, where now the squeezing (stretching) effect is promptly verified in the 
discrete variable v {p). In both cases, the appearance of the squeezing and stretching effects can be understood in this 
context as a direct consequence of the construction process adopted for the discrete representatives of the states under 
investigation. Indeed, the introduction of an extra parameter in the second argument of the Jacobi thcta functions has 
allowed us to properly construct out, from the finite-dimensional context, a wave function that inherently embodies 
the discrete analogues of the desired properties of the continuous squeezed states. 

However, it is worth noticing that we have only discussed hitherto the action of a hypothetical squeezing generator 
on the discrete vacuum state without deriving this operator and/or establishing its mathematical properties. This 
breach will be filled in the appendix B through a discussion on the discrete representative of the squeezing generator 
for finite Hilbert spaces. 



B. Discrete representatives of the squeezed coherent states 



Initially, let us consider the decomposition of the diagonal projector related to the discrete coherent states 
{\p.,v)}ji,u=-t,...,i, that is 
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Since the discrete representatives of the squeezed coherent states are defined within this context through the action 
of the squeezing operator 3C(s) on the discrete coherent states, the diagonal projector 

e 

\fi,D; S )(fi,P; S \= Tt[sV>Olft*;*><fcM]S(T?',0 (14) 

n',Z'=-i 

can be interpreted as being the decomposition of 3t(s) |/2, v}{ft, v\3t'(s) in a complete orthonormal operator basis 
characterized by a set of iV 2 -operators {S(r/, £')}, where the coefficients (p,, v\ s\S* (rj' , £')\p,, v\ s) are calculated through 
the formal expression 

JV-l 

N r „ 1 (-i>,^;l)r„ m (-r; , ,-C';s)r„o(^,-Ai;l) (15) 

m,n— o 

with the nondiagonal matrix elements in the finite squeezed-number basis T mn {a, f3;s) '■— (m;s\S(a, /3)\n;s) defined 
as Q3 

r m „(a,/3;s) = exp ^ exp ^^. K aj C^,mW^,nW- 

It is worth mentioning that these generalized nondiagonal matrix elements satisfy the relations r mn (0, 0;s) = 5m} n 
and roo(a,/3;s) = J^(a,/3), which are related to the orthogonality rule for the finite squeezed-number states and 
the diagonal matrix element (0;s|T( s )(^i, ^)|0;s) for the squeezed-vacuum state. Furthermore, for m = n and 5 = 1 
it is possible to verify that T nn (a, j3; 1) coincides with J(f n (a,(3). Next, we will derive an expression for the overlap 
probability distribution P„(/l,P;s) by means of an alternative form to that established in this paragraph for the 
diagonal projector (fT4|) . 
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The overlap probability distribution P n (/i, v\s) := \(n\fi, v;s)\ 2 associated with the discrete representatives of the 
squeezed coherent states can be evaluated directly from Eq. (TB|) as follows: 
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On the other hand, if one considers the mathematical procedure used to derive Eq. (jTTJ) in the context here exposed, 
the alternative expression 
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allows us to show that P n (p,9;s) can be also estimated through a product of the kernel-like function K(p, 9, r/, s) 
- which describes how the squeezing effect occurs in an A^-dimensional phase space labeled by the set of discrete 
variables {p, 9} - and the double Fourier transform of the initial distribution function P n (p! , z/), 



P n (p, v;s) 
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exp 
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K(H,9,r)',e;s)P n (ii',v'). 



In the absence of squeezing effect, K (/2, 77', 1) coincides exactly with exp [— (2ni/N)(r]' fi + and consequently, 
the double sum over 77' and f leads us to reobtain, by means of the product of Kronecker delta functions S^^S^p, 
the discrete probability distribution function ©. 



IV. APPLICATIONS 



In Ref. [3l| it has been mentioned the concrete possibility of using the degenerate parametric amplification of a 
signal to reproduce, from the experimental point of view, the unitary transformation <S^(£)p<S(£) gene rated by the 
action of the squeezing operator upon physical systems described by continuous quantum variables [32j. However, if 
one deals with specific physical systems characterized by quantum states of arbitrary finite dimension and described 
by discrete quantum variables, the squeezing action will be modelled hypothetically through a nonlinear medium or 
other external apparatus connected to the system of interest. This connection will then produce a global effect upon 
the system that consists in modifying its discrete wavefunction according to the rules established by a particular 
squeezing operator. In this sense, let us suppose that 3C(s) can be modelled by one of these hypotheses which 
leads one to produce, within certain experimental and theoretical limitations (e.g., the degrading and ubiquitous 
decoherence due to the unavoidable coupling with the environment), the unitary transformation 3£^(s)pj£(s) for 
physical systems belonging to finite-dimensional Hilbert spaces. Next, we will derive an analytical expression for 
F^(p,9;s) ■= Tr[T( s ) (p., 9)3£^ (s)p3£(s)} which allows us to obtain, as particular situations of this parametrized 
function, the discrete Husimi (s = —1), Wigner (s = 0), and Glauber-Sudarshan (s = 1) functions in an iV 2 - 
dimensional phase space. 

For this purpose, let us initially consider the cyclic invariance property under trace operation in the definition of 
F( s )(/2, 9;s), namely 

F^(p,,9;s) = Tr[X(s)T^(p,9)X' i (s)p}. 

The next step consists in substituting Eq. (JTJ) on the right-hand side of this equality. Thus, after some calculations, 
we immediately recognize 



F^(p,9; 5 ) = ^= exp 



s (s) (?7,£;s) 



(18) 



as being the desired analytical expression for the parametrized function where, in particular, 

sW(r 7) C;5) = Tr[X( S )SW(r ? ,03e t ( 5 )p] 
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FIG. 2: Generalized version of the scattering circuit used to evaluate the real and imaginary parts of the expectation value 
Tr[U3£ t (s)/3j£(s)] for a unitary operator U = \^NS(ri, £), where |0) represents the ancillary qubit state which acts as a probe 
particle in a scattering experiment, and H denotes a Hadamard transform. The measurements of the ancillary qubit polariza- 
tions along the z and y axes allow us to construct the discrete Wigner function W(/i,i/;s) = Tr[T*°'(/i, v)^ (s)pX(s)] (discrete 
characteristic function H (0) (r?,£;s) = Tr[S (0) (?7, ^)3£ t (s)p3£(s)]) in the presence (absence) of the controlled-FT operation. 



represents the squeezed discrete s-ordered characteristic function. Note that for s — 1, the parametrized function 
F^(jl,D; 1) coincides exactly with that obtained by Eq. Now, we will concentrate our efforts in describing a 

physical process which allows us to measure directly the discrete Wigner function W(/2, D;s) for any quantum system 
constituted by a finite-dimensional state space. 



A. Measuring the discrete Wigner function via generalized scattering circuit 

Recently, we have employed a slightly modified version of the scattering circuit proposed by Paz and co-workers (23j 
to measure the discrete Wigner function W(/i, v) in the absence of squeezing effects [18j . In the present approach, we 
generalize such version in order to introduce within the primary circuit an auxiliary set of programmable gate arrays 
which allows us to reproduce the unitary transformation 3& (s)pX(s). This theoretical implementation will then open 
new possibilities of measurements on the generalized circuit: for instance, both the discrete Wigner function W(/2, P; s) 
and the characteristic function E^ -* (77, £; s) will be inferred in the same experiment through a simple operation associted 
with the controlled Fourier transform (FT). 

Following, let us describe some important implementations introduced within the scattering circuit which is il- 
lustrated in Fig. 2. Basically, we modify this circuit by inserting a controlled-U operation between the Hadamard 
gates, with U = VWSfa, £) acting on a squeezed quantum system whose density operator is exactly described by 
X\s)pX(s), and also a controlled-FT operation after the second Hadamard gate. This procedure allows us to obtain 
a set of measurements on the polarizations along the z and y axes of the ancillary qubit |0) that leads to the expec- 
tation values (cr 2 ) = y/NKe\W(p,,D;s)] and (er y ) — \A/VIm[W(/I, v;s)], respectively. Furthermore, in the absence of 
the controlled-FT operation, these measurements permit us to reach the squeezed characteristic function fa, £; s) 
for s = - i.e., (cr z ) = v^Re^f^^s)] and (cr y ) = v^VIm[S(°) fa, £;*)]. 



B. A theoretical study on the quantum interference effects 



To conclude this section, we will present a theoretical study on the quantum interference effects between the discrete 
variables of an 7V 2 -dimensional phase space. Indeed, the concept of information entropy for the continuous squeezed 
states [33[ will be extended in order to include a consistent mathematical formulation based on the discrete squeezed 
states. Hence, let us introduce the entropy functional 

E[W;s]:=-l £ W(/i,i/;s)In[W(/i,i/;«)] , (19) 

H,v=-i 

where H.(p, v\ s) ■= (fi, v\ s|p|/J, v\ s) represents the discrete Husimi function related to the displaced squeezed- vacuum 
states. The virtues of this definition are directly associated with its continuous counterpart, that is, (i) once 7i(/x, v\5) 
is a positive definite probability distribution function and limited to the interval [0,1], the joint entropy E[7i;s] 
characterizes a well-defined function whose behaviour does not present any mathematical inconsistencies. Besides, 
(ii) this joint entropy essentially measures the functional correlation between the discrete eigenvalues of coordinate 
and momentum-like operators for a specific class of quantum systems with finite Hilbert spaces. Finally, it is worth 
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mentioning that E[W;s] can also be defined for the discrete representatives of the squeezed coherent states or any 
discrete representation related to a finite phase-space. 

Next, let us consider the functionals related to the marginal entropies 



E[Q;s] 



1 



]T Q(/i;s)ln[Q(/i;s) 



(20) 



and 



1 1 



(21) 



which depend on the discrete marginal distributions [if 



27ri \ ^foO) p( _ x) 



27ri \ jr g (o,fl „ ( _ 1} , 
iv 5 ; jt(o,o 



^; S ) = ^exp(-^]^^Sl-"(0.O 



Consequently, the conditional entropies are given by 



E[H/Q;s}:=-- ]T H{h,v\b)\b. 



n,v=-e 



and 



E[H/K;s] 



N 



E H(p,v\a)\n. 



'H{n,v\s) 

C(w«) 



= E[H;s]-E[Q;s] 



E[H;s] -E[K;s]. 



(22) 



(23) 



Note that (j22|) and (I23[) allow us not only to establish the balance equation 

E[H/Q; 5} + E[Q; s] = E[H/U\ s] + E[TZ; s] , (24) 
but also to derive some mathematical relations among these functionals from the Araki-Lieb inequality [3~ij 



\E[Q;s]-E[K;s]\ < E[H;s] < E[Q;s] +E[K;s]. 



(25) 



For instance, it is easy to show that the inequalities E[H/Q;s] < E[lZ;s] and E[H/lZ;s] < E[Q;s] can be promptly 
obtained from the definitions established for the conditional entropies and the right-hand side of Eq. (|25|) - which 
defines the subadditivity property. Moreover, the equal signs hold in both situations only when the discrete variables 
{/Lt, u} e [— £, £] are functionally uncorrelated, namely, v,s) factorizes in the product Q(/i; s)H{v; s) . 

To avoid any ambiguity in the significance of the subadditivity property and to clarify its precise meaning, it is 
useful at this moment to define the correlation functional 



C[H;s] 



E[Q;s] +EpK;s] -E[H;s] 
E[Q;s] +E[^;5] - Min(E[W; s}) ' 



(26) 



where Min(E[7i ; s] ) denotes the global minimum of E[Tt; s] for a given dimension of the Hilbert space and any squeezing 
parameter [35|. Essentially, C[W;s] measures the strengh of the functional correlation between the discrete variables 
fi and v, having as reference a factorizable Husimi function 7i(/i, j/;s). In addition, Eq. (|26[) assumes any values 
restricted to the interval [0, 1] such that its inferior limit corresponds to the situation where the joint distribution 
function factorizes, while the superior limit is associated with the opposite situation - i.e., when the discrete variables 
are completely correlated. Hence, C[H;s] can be considered as a first reasonable measure of functional correlation 
between the discrete variables of an A^ 2 -dimensional phase space, being its extension to multipartite systems very 
useful within the context of quantum information theory [24| and/or other possible scenarios in physics |36| . 

As a first practical application for these entropy functionals let us consider the discrete Husimi distribution function 
Pq([a,v]s) associated with the lowest level of the finite number state (namely, the vacuum state). Figure 3 shows 
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FIG. 3: Plots of some entropy functionals associated with the lowest level of the finite number state (vacuum state) as function 
of the dimensionless squeezing parameter s for N = 3 fixed. Picture (a) shows the behaviour of the joint entropy E[Po;s] 
related to the discrete Husimi distribution function Po(fi, v; s) , where its minimum value for s = 1 is approximately given by 
0.625948. The solid and dashed lines represent in picture (b) a graphical illustration of the marginal entropies E[Qo;s] and 
E[72.o;s], respectively. Note that both marginal entropies present a distinct behaviour when s — > or s 2> 1, and this fact 
implies in the absence of knowledge about one of the discrete variables /i or v. The correlation functional is then plotted in 
picture (c) and its maximum value, attained in a = 1, corroborates the numerical results obtained in the previous figures. 



the plots of (a) joint entropy, (b) marginal entropies, and (c) correlation as function of the squeezing parameter for 
N = 3 fixed. In particular, we verify through numerical investigations that the minimum value attained by E[Po;s] 
in picture (a) for s = 1 depends essentially on the Hilbert-space dimensionality. To illustrate this kind of dependence, 
we have evaluated the joint entropy E[Po; 1] for N = 3, 5, 7, 9 and obtained, as final results, the approximated values 

0. 625948, 0.953965, 0.992272, and 0.998598, respectively. Furthermore, the asymptotic value Min(E[P ; 1]) -> 1 may 
be approximately reached for N > 11, this value being coincident with that found for the continuous analogue (33[. On 
the other hand, it is worth noticing that E[Qo;s] (solid line) and E[72.oi-s] (dashed line) in picture (b) have a different 
behaviour when s -> or s > 1, and this fact can be explained by means of the absence of knowledge (uncertainty) 
about the discrete coordinate and momentum- like variables for each specific situation of s. Since the minimum value 
of E[Po;-s] implies in maximum information gain on both the discrete variables of a finite-dimensional phase space, it 
is expected that C[Pq;s] attains a maximum value in this case (indeed, this is what exactly happens in picture (c) for 
5 = 1). However, if one considers the continuum limit A 7 ^ — * oo, it is easy to show that Po(q,p] A) = Qo(q] A)i?o(p> A), 

1. e., the Husimi distribution function is completely uncorrelated, and consequently, the correlation functional becomes 
zero for any value of squeezing parameter. 



V. CONCLUSIONS 



In this paper we have established a set of important formal results that allows us to obtain a discrete analog for 
two well-known definitions of squeezed states in a continuous phase-space (37l [38| . In fact, we have employed the 
mathematical fundamentals developed in Ref. [l8| with the aim of constructing out a consistent formalism for the 
discrete representatives of the displaced squeezed- vacuum states and squeezed coherent states through their diagonal 
projectors in an iV 2 -dimcnsional phase space. Next, we have applied our formalism to the context of quantum 
tomography and quantum information theory in order to obtain a finite phase-space description of an important 
topic in physics which was nominated in the recent past as quantum interference effect (39j . It is worth noticing that 
such description has allowed us to attain new results, within which some of them deserve to be mentioned: (i) we 
have introduced a generalized version (where the squeezing effects are taking into account) of the scattering circuit 
proposed by Paz and co-workers [23] to measure the discrete Wigner function or its associated characteristic function 
for any physical system described by an iV-dimensional space of states; (ii) we have also established a first reasonable 
measure of functional correlation between the discrete variables of a finite phase-space whose extension to quantum 
information theory, as a measure of entanglement in multipartite systems, is completely feasible; and finally, (iii) we 
have performed numerical investigations involving the discrete squeezed- vacuum state and showed that, in particular, 
the maximum value of this measure presents a strong dependence on the Hilbert-space dimension. As a concluding 
remark, it is worth mentioning that the formalism developed here opens new possibilities of investigation mainly in 
modern research on quantum optics as well as the foundations of quantum mechanics. These considerations are under 
current research and will be published elsewhere. 
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APPENDIX A: CALCULATIONAL DETAILS OF THE DIAGONAL MATRIX ELEMENT (n|S(?7, £)\n) 



In this appendix, we will explore the technique of breaking infinite sums in mod(TV) equivalence classes to derive a 
closed-form analytical expression for the diagonal matrix element J^ n (r],^) = y/N (n\S(r},£)\n), with 



•%n(v,Q = ex P ( -Tfyt) exp ( ~W ari ) 



(Al) 



written in terms of the coefficients [ll|, [26 



N, 



(-i)« 



N„ 



^(7ro)"iV 93' 



■ exp(-7m3 2 ) i?3 [2o(«; — 3) |2ia] 



(A2) 



3=0 



where N n is the normalization constant and H n {z) is a Hermite polynomial. Basically, our goal will be split up into 
two parts: the first one corresponds to the evaluation of the normalization constant N„, while the second one is 
associated with the derivation of the generalized function J^n,(rj, £). 

1. Normalization constant 

The constant N„ can be evaluated directly from the normalization relation for the finite number basis 
{\n)}„=o,...,N-i, namely 



1 = (n\n) = ^ {n\u K )(u K \n) = ^ K,n$K,- 



(A3) 



K=-l 



where {\u k }}k=— represent the eigenstates of the Schwinger unitary operator U. Then, substituting the expression 
for n in Eq. (|A3|) we obtain, as a first step in the calculations, the equality 



N 2 

_ ] n 

N 



J2 J2 exp 



K =-£ {a,/3}& 



N 



N 



H n I Jjf<*\ //„ ( \l=±-l 



At this point, it is worth noticing that the sum over k can be readily carried out, which gives us NS^l (the superscript 
[TV] on the Kronecker delta denotes that this function is different from zero when its labels are congruent modulo N). 
Thus, a will assume the values (3 + (3'N, with arbitrary 0' G Z, yielding 

1 = N 2 £ cxp {-^[,3 2 + (p + pN) 2 ]}H n (Mp)H n \M(P + / 3>N) 

{/3,/3'}GZ ' \ V / L V 

The second step consists in considering the generating function for the product of two Hermite polynomials [4C 



H n (x)H n (y) = — (l-4 3 2 )- 1 / 2 exp 



4xj/3 - 4(.t 2 + y 2 )f 



3=0 



in the last equality, which permits us to reach, after some manipulations, the intermediate result 



i = n 2 ^ r (i-4 3 2 r i / 2j(3) 



(A4) 



3=0 
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with J(j) given by 



J(3) = ^ exp 

{/3,/3'}ez 



2tt 1 + 4 3 2 (o FN 
N 1 - A^ 2 \ 2 



ttN 1 + 4 3 2 , 2 2tt 4 3 



2 1-43' 



N 1 — 43 



;{3(f3 + f3'N) 



Note that the sum over /3' can be separated out in two contributions coming from the even (e) and odd (o) integers, 
and consequently this procedure implies that J (3) = J e (3) + J D (3). Now, let us go one step further in order to determine 
each term separately. 

For example, the even term J e (3) can be dealt with by shifting the sum over (3 by —/3'N. This fact produces a 
decoupling between the discrete indices /3 and (3', namely 



Je(3) = E exp 



2?r 1 - 23 



HE 



exp 



which allows us to identify each sum with a particular Jacobi theta function as follows [251 ]: 

Je(3) = H/(3)] 1/2 ^[0|4/(3)]^ 3 [0|/(3)], 



(A5) 



where f(i) — ia(l — 2g) (1 + 23) 1 and a = (2N) 1 . The odd term J Q (3) can also be determined by means of a similar 
mathematical procedure, yielding as result the analytical expression 



J o (3)-hi/(3)] 1/2 ^[0|4/(3M[0|/(3)]. 



(A6) 



Therefore, with the help of (|A4[) . the sum of these contributions leads us to obtain a closed form for the normalization 
constant N„. 



2. Generalized function 



Initially, let us substitute the analytical expressions of the coefficients $cr,n and 5"* _^ n into Eq. (|Al|) . The sum 

over a yields the result NS^0 + , which implies that a will assume the values (3 + rj + (3'N for (3' £ Z. Thus, with the 
help of the generating function for the product of two Hermite polynomials and after some algebraic manipulations, 
the generalized function can be written as N„.^ n (r/, £), where ^n(V->£) is defined through the equation 



.fa>0 gpra 



2 3 )- 1 M(/ ? ,e;3: 



3=0 



(A7) 



and M(?7,£;3) denotes an auxiliary function given by 



M(ry,£;3) = [-1/(3)] 



-1/2 



exp 



x ex P 



2tt 1 + 43 2 
' N 1 - 43 2 



N 



7T 1 



AT I-43 



2 

2»r 



4 3 9 



ttN 1 



2 I-43 



E ex p 

43 V' 2 4 



2tt 1 + 43 2 
N 1 - 43 2 

2^ 43 



77/3- 



2tti 
~N 



N l-Af 



((3-r))((3 + (3'N) 



(A8) 



The next step consists basically in evaluating the sums over (3 and (3' by means of the mathematical procedure adopted 
for J(3) - i.e., the auxiliary function will be decomposed in two parts coming from the even and odd contributions 
associated with the discrete index (3'. 

It is worth mentioning that the technical details involved in the determination of M(?7, £; 3) = M e (ij, £; 3) +M (t], £; 3) 
will be suppressed at this moment without significant consequences for our initial purpose, being exhibited only the 
most important results. In this sense, the even contribution 



Me(»?,63) = H/( 3 )] 



-1/2 



exp 



i7T 7T 1 



2 



43 _ a 



exp ( -2itN- 

/3'ez ^ 



V 



23 



2tt 



43 

I + 23 
1-23 



2^ exp 

/3GZ 



2tt 1 



AT I + 23 



V 



77/3' 



27ri 



I-23 
1 + 23' 



13 



assumes the simplified form 

while the odd contribution 

M (77,£;3) 



{^3[a/?|/(3)]^3[<|/(3)] + t9 3 [ar7|/(3)]^4[aC|/(3)]e i7rr '} , 



x J^exp 

/3ez 



2tt 1 



1/2 



exp 
23 



i7T 7T 1 


AT /S TV 1 - 4j 



2 

2 T 



4 3' _2 



l-2 3 

' 1 

1 + 23 



E 

/3'GZ 



exp 



iV 1 + 23 

l + 2 3 



N 
~2 



2 



2?ri 



~2tt7V- 



I-23 



N 



2tt 



,1-23 

L 7 

1 + 23 



77/? 



I + 23 
j 

I-23 



can be represented through the expression 



M (^;3) = \ {Mav\f(i)¥sH\f(i)]e^ + Mav\f(^M\f(i)}e^ +N) } . 

Note that some fundamental properties of the Jacobi theta functions were extensively used in the derivation process 
of M e ( D ) (77, £; 3), such properties being promptly listed and discussed in Ref. [25j. In addition, the sum of these 
contributions leads us not only to determine the function ^ n (rj, £) through the mathematical operation defined in Eq. 
(|A7|) . but also to show that the normalization constant N n is connected with ^#„(0, 0) by means of the equality N~ 2 = 
.«# n (0, 0). Consequently, the generalized function can be written in this context as J^niViO = -^n ( T L 1 -^n (0,0), 
which allows us to obtain directly the diagonal matrix element (n|S(?y, £)|n). 



APPENDIX B: THE SQUEEZING GENERATOR FOR FINITE-DIMENSIONAL HILBERT SPACES 



We initiate this appendix establishing some particular results for the squeezing generator S{r) in infinite-dimensional 
Hilbert spaces which constitute the continuous counterparts of the discrete cases to be presented in this section. For 
instance, the normalized wavefunction for the squeezed coherent states can be obtained through a simple mathematical 
procedure sketched out by the relation 



x\q,p;r) = / dx'{x\S(r)\x'){x'\q,p), 



where 

(x'\q,p) — n "' ' exp 

represents the wave function for the continuous coherent states in the coordinate representation and 



— 7T ' ' 



V'[X - - 



\S(r)\x') = A 1/4 exp 



S(x' - X^x) 



(Bl) 



(B2) 



(B3) 



is the nondiagonal matrix element associated with the squeezing operator S(r) when A := cxp(— 2r). The adequate 
manipulation of these results allows us to show that (jBip can then be promptly derived by means of the identity 
(x\q,p;X) — X 1 ^ 4 {X 1 ^ 2 x\q,p). In fact, this simple operation leads us to obtain a generalized set of wave functions 
belonging to a wide class of squeezed states of the electromagnetic field characterized by [H 

{|¥;r) :=S(r)|*) : 5(r) = exp[-i(r/2){Q, P}], re (-00,00)}, 

namely (x|*;A) = A 1 / 4 (A 1 / 2 X |*). 

Now, let us apply this result to a particular set of nonclassical states described by the number states {|n)} n pN- 
Thus, the normalized wavefunction associated with the squeezed number states assumes the analytical form [41| 



(x\n;X) 



A 



7r(2"n!) 2 



1/4 



1 



exp[-^)H n (X 1 ^x), 



(B4) 
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whose complementary wavefunction in the momentum space representation, 

^' n;A > = ; w^ 1 !" a/4 ex P (~\ x ^pi) H n (\- 1/2 Px ), (B5) 

[7rA(2™n!) 2 ] 1 V 2 / 

permits us to establish a well-known important link between both the coordinate and momentum representations via 
the Fourier transform, 

(x\n;X)=f -^exp(ixp x )(p x |n;A). (B6) 



Basically, these results will be our guide-lines in the construction process of the discrete counterpart for the squeezing 
generator in finite-dimensional Hilbert spaces. 

In appendix A we have shown that are the coefficients that lead us to connect the finite number states 

{\n}} n =o,...,N-i and the discrete coordinate-like representation {|u K )} K= _£,.../ through the expansion 

t 

|n) = Sk,™K). (B7) 

K=-l 

Following, let us define the squeezed number states in the context of finite-dimensional Hilbert spaces by means of 
the mathematical statement 

t 

\n;s) := X{s)\n) = ^ S K ,„(s)K), (B8) 

K=—l 

where {^ K , n (s)} represent the generalized coefficients expressed as 



y/(iras 2 ) n N 3j 
N„(s) being the normalization constant given by 



Xia) 9 " exp(^a5V)^ 3 [2a( K -3)|2ia S - 2 ] 



(B9) 



3=0 



N- 2 (s) = ^_(l + 2 3 )- 1 {tf 3 [0|4 S - 2 /(3)]^[0| S 2 /(3)] +^[0|4 6 - 2 /( 3 )] l ? 4 [0| S 2 /(3)]} 



3=0 



Note that the squeezing parameter s plays, in principle, the same functional role reserved to the continuous counterpart 
A - compare Eqs. (|B4[) and (|B9p - i.e., both the squeezing parameters modify the widths of their respective wave 
packets and this fact allows us to introduce satisfactorily the squeezing effects within the context of finite-dimensional 
Hilbert spaces. In addition, the wavepacket for the discrete squeezed-vacuum state © can be considered as a particular 
component of the generalized coefficients since it corresponds to the specific case 3k,q(s). Consequently, the complete 
set of components characterized by {3 Kl n(s)} n= o,...../v_i will be sufficient in this approach to determine precisely the 
squeezing operator 3E(s) and its inherent mathematical properties. 

The action of X(s) on the finite number states not only defines the squeezed number states (|B8|) but also leads us 
to establish a formal expression for the squeezing generator as follows: 

N-l 



X(s) := IMH- (BIO) 



In what concerns the operator X(s), it is worth noticing that its representation in the set of Mehta's states [26| carries 
a peculiar behaviour inherent to those states. Indeed the study of this particular representation deserves a detalied 
investigation since the peculiarities mentioned before come from nonorthogonality of the Mehta's set - for technical 
details on theoretical and numerical investigations of this specific topic, see Ref. [27]. Finally, let us mention briefly 
that the proof on the unitarity of the discrete squeezing operator will be published elsewhere [42j ■ 
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